Linear response to perturbation of non-exponential renewal processes 
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We study the linear response of a two-state stochastic process, obeying renewal condition, by 
means of a stochastic rate equation equivalent to a master equation with infinite memory. We show 
that the condition of perennial aging makes the response to coherent perturbation vanish in the 
long-time limit. 
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Many complex processes generate erratic jumps back 
and forth from a state "on" to a state "off'. We limit 
ourselves to quoting ionic channel fluctuations 0, 0, 0] , 
currently triggering the search for a form of stochastic 
resonance valid also in the non-exponential case 0, and 
the intermittency of blinking nanocrystals 0, @ : It has 
been assessed that the intermittent fluorescence of these 
materials obey renewal theory 0, namely, a jump from 
one to the other state, has the effect of resetting to zero 
the system's memory. The non-exponential nature of the 
distribution of sojourn times makes this renewal process 
non ergodic and generates aging effects that are the ob- 
ject of an increasing theoretical interest 0EJ. Similar 
properties are found with surface-enhanced Raman spec- 
tra of single molecules |10j. 

The authors of proved that aging, in conflict with 
the pioneering work of Ref. |l2j | claiming the equivalence 
between Generalized Master Equation (GME) and Con- 
tinuous Time Random Walk (CTRW) [l|, makes the 
GME response to perturbation different from the CTRW 
response. These authors, however, did not propose a new 
approach to the linear response. 

Here we present an approach to the linear response 
valid in the aging as well as in the stationary case, as in 
earlier work [j, [lj] . The key ingredient of the method is 
to express the CTRW prescription through a Stochastic 
Master Equation (SME), which yields the linear response 
under the natural assumption that the time distribution 
of collisions is not significantly affected by perturbation. 

Let us consider a two-state renewal process, and, for 
simplicity, let us assume that the distribution of times of 
sojourn in the state "on" is the same as the distribution 
in the state "off" . We assign to the Survival Probability 
(SP) of this process, ^f(t), the inverse power law form 



with ji > 1. This corresponds to the joint action of the 
time dependent rate 0, r(t) = r /(l + rit), with 
ro = (/i — 1)/T and r\ = 1/T, and of a resetting pre- 
scription. To illustrate this condition, let us imagine the 



random drawing of a number from the interval / = [0, 1] 
at discrete times i — 0,1,2.... The interval / is divided 
into two parts, I± and I2, with I\ ranging from to pi, 
and I2 from pi to 1. Note that pi = 1 — qi < 1 and 
qi << 1, and, as a consequence, the number of times we 
keep drawing numbers from I\ , without moving to I2 , is 
very large. 

Let us evaluate the distribution of these persistence 
times, and let us discuss under which conditions we get 
the SP of Eq. (|TJ). The SP function is the probability 
of remaining in I\ after n drawings, and is consequently 
given by 
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Using the condition qi << 1, and evaluating the loga- 
rithm of both terms of Eq. J5J), we obtain 

n 

M*(")) = -£ft- (3) 

i=l 

The condition q t « 1 implies that i and n of Eq. @ 
are so large as to make qi virtually identical to a function 
of the continuous time t, q t = q(t) = rorjit), with rj(t) = 
1/(1 + nt). Thus, Eq. © yields the SP of Eq. ©, 
and the corresponding waiting time distribution density, 
ip(r), reads 

rpn-l 

^) = (M- 1 ) F TTF- W 

We denote as collisions the rare drawings of a number 
from I2, followed by resetting. Thus the collisions occur- 
ring at times t\, ti+t 2 , yield: rj{t) = 1/(1 + rqt), < 
t < n; rj(t) = 1/(1 + rx(t - n)),rx < t < n + r 2 , and so 
on. Note that ?/(0) = 1 means that we prepare the sys- 
tem at time t = 0. We adopt a coin tossing prescription 
to decide whether to keep or to change sign, after any 
collision. 

The state of the system after the n-th collision is de- 
scribed by the two-dimensional vector P(n), whose com- 
ponents, P\ and P2, are the probabilities of finding the 



2 



system in the corresponding states, 
produces the change 



The n-th collision 



P(n) = MP(ii - 1), 



where 



M 



1/2 1/2 
1/2 1/2 



(5) 
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Let us denote by ui(i) and 02 (t), the probabilities of 
finding the system in the state "on" and "off', respec- 
tively, at a time t, which, does not necessarily correspond 
to the collision occurrence. These probabilities are driven 
by the following SME 



d a 2 (t) 
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-^i(*) 



^i(t), 



(7) 



with r+(t) = r_(f) = r(t) = r f(f). The condition 
r+(t) 7^ r -(£) is hereby used to discuss the response to 
perturbation. The variable £(t) always vanish but in the 
correspondence of a collision, where it gets the value 1 /ro, 
so as to make Eq. Q equivalent to Eq. (J5J, Note that 
time t is discrete and that the derivatives dai/dt, with 
i = 1,2, are defined by dai/dt = Of(i + 1) — <7j(t). In con- 
clusion, the collision time is a stochastic variable making 
the rate r(£) fluctuate between and 1. 

To derive the linear response to perturbation, it is nec- 
essary to connect the SME of Eq. to the CTRW 
functions *&(t) and tf)(t), including the corresponding ag- 
ing properties as well. To establish this connection, let 
us notice that the CTRW picture can be converted into 
the non-Markov GME [3 (13 



jP(t) = - J dr$(t - r)KPp(r) 



The matrix K is defined by 



K = 



1 



(8) 



(9) 



and the memory kernel is related to tp(t) in the 

Laplace domain through 



uip(u) 
l-$(u)' 



(10) 



where $(w) and ip(u) denote the Laplace transforms of 
<!>(£) and ip(i), respectively. 

Now, let us define the quantity 



n(t) = Px{t) - p 2 (t) = 2Pi(t) - 1. 



(ii) 



Using Eq. (JSJ, we get for the Laplace transform of LI(£) , 
n(it), the following expression: 



n(u) = 



(12) 



which, in turn, using Eq. 
namely, the key property 



(fTOf) . gives fl(u) = *(u), 



IL(t) = (t). 



(13) 



Let us go back to Eq. Q with no perturbation. The 
average on the fluctuating rate yields < <Xj(£) >= Pi(t), 
with i = 1,2, and the average of the stochastic quantity 
E(t) = <J\{t) — <J2(t) = 2<7i(i) — 1 becomes equal to II(i) 
of Eq. llljl . To double check this important prediction, 
let us assume that all the systems are located at the be- 
ginning of their sojourn in cither the state "on" or the 
state "off" , so as to fit the prescription of preparing the 
system at t — 0. Thanks to Eq. (JJJ, we obtain 



£(£) = exp 



-To 



at')df 



(14) 



which yields 



(£(*)> = *(*)• 



(15) 



In fact, E(£) changes only at the occurrence of a collision, 
and \&(£), as we have seen, is the probability that no 
collision occurs up to time t. Using Eq. (T5jl. we conclude 
that < E(t) >= U(t). 

With the help of Eq. Q we establish the condition 
for the linear response to occur. Let us assume r±(t) = 
ro£(i)(l + eF±(t)), where e is the perturbation strength 
and the functions F±(t) describe the action of external 
perturbation on the corresponding states. Using Eq. Q 
we obtain 



dt 



m = -r„£(t)(l + eS(t))S(t) - r S(t)ef(t)), (16) 



where S(t) = (F + (t) + F_(t))/2 and f(t) = (F + (t) - 
F_(£))/2. The requirement that the system response 
does not vanish enforces the condition f(t) ^ 0, and 
the additional request that this response be linear yields 
eS(t) << 1. Note that the earlier illustrated ideal numer- 
ical experiment, in the presence of perturbation, should 
be done with q±(t) = r r)(t)(l + eF±(t)). With q±(t) re- 
maining very small, the occurrence of long-time collisions 
is determined by very small values of T)(t), with a very 
small dependence on F±(t), which is sufficient, however, 
to produce a time-dependent bias on LI(£) . This makes 
it natural to assume that £(£) is independent of pertur- 
bation. 

Using the linear response approximation and solving 
the corresponding approximated form of Eq. i|16|) , we ar- 
rive at 

E(t) = -r Q e J" dt' exp (-r„ jf £(*")<**") £(*W) 



^(0)exp [ -r J Z(t')dt') 



(17) 
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We are interested in the mean value of £(t), namely II(i). 
This leads us to write 



n(t) = -e / x(t, t')f(t')dt' + n(o)*(t), (is) 

Jo 



where 



(19) 



with ^(f , t') being the characteristic function 

s ^(-r jf dt"£(t"))) • (20) 

With 11(0) = 0, Eq. l(T^|l becomes formally identical to 
the common linear response prescription 18], with ag- 
ing violating, however, the stationary condition %(t, t') — 
X (t-t>). 

Using the same arguments as those earlier adopted to 
prove Eq. I|15|). we show that \I/(i, t') coincides with the 
survival probability of age t' of Ref. 19] . According to 
[l9| . the function ty(t,t') coincides with the aging cor- 
relation function of the dichotomous signal under study, 
C(t,t'). By means of the stationary assumption we re- 
cover the results of Refs. 0, and we realize why 
these earlier results go beyond the Green-Kubo predic- 
tion mini. 

In the case /i < 2, the stationary condition cannot be 
realized, not even in principle, while in the case /i > 2 it is 
possible, even if /i < 3 makes very slow the relaxation to 
equilibrium p^ . Let us study first the more traditional 
condition /i > 2, and let us switch perturbation on after 
the transient process necessary to reach equilibrium. We 
have El 



x(t,t') = -i/j 00 {t-t') = -(fx-2) 



rpjj,- 



(t-t' + ry 



- (21) 



the stationary correlation function of the "on" and "off' 
fluctuation, C(t — if), exists, with the analytical form 
(T/(f + T)y- 2 , thereby yielding 



i> as {t-t') = --C{t-t'). 



(22) 



Note that the linear response of Eq. (fT%)l . with — x(t, t') 
replaced by ipoc(t — t') of Eq. (|22|l coincides with the 
prescription of the phenomenological approach of Ref. 
[4J. In the case where f(t) = @(t), with 0(i) denoting 
the unit step function, by plugging Eq. (|2"2")l into Eq. 
(jT%)l . we obtain 



n(oo) = e, 



(23) 



and consequently the constant drift of Ref. [lj] . 

In the case where f(t) = Q(t) cos(uit), namely the case 
of stochastic resonance discussed by the authors of Ref. 
, by using the expression for the Laplace transform of 



tpooit) and the convolution theorem, we obtain for the 
Laplace transform of II (t) the following result: 



H(u) = e(l + cT^- 2 ^" 2 ) 



(a 2 



(24) 



where c = — L(3 — fi). Using the fractional derivative 
method of Ref. [24j , we evaluate the anti-Laplace trans- 
form of tl(u) of Eq. l|2"3jl. In the case ujT « 1, we 
obtain [2^, |23|: 



n(t) 



1 + (uTf 2 T (3 - n) cos -fj, 



cos cut 



-e(u)Ty 2 r (3 — n) sin — /i sinwt. 



(25) 



As earlier pointed out, the theory of this Letter applies 
also to the aging case, and to the condition \i < 2, which 
makes aging a perennial condition of renewal systems. It 
is known [23 (see also Ref. 0) that the exact expression 
for the aged tfj(t,t') is 



00 „t 

V(m') = <K<) + ]T / 

n=l Jo 



ip„{T)ip(t - r)dr, 



(26) 



where ^> n {t) denotes the probability density of occur- 
rence, at time t, of the last of a sequence of n collisions. 
Let us go back to Eq. (|T%)| and let us set 11(0) = 0. 
This has the effect of producing a condition where a non- 
vanishing II(t) appears only as a perturbation effect. For 
the Laplace transform of LT(t), we obtain the following 
expression 



n(u) = -eRe(E(u)), 
where E(u) is the Laplace transform of 



E(t) 



dt't(j(t,t')exp(~iujt'). 



(27) 



(28) 



After some algebra, we find 



• \ i}){u + iuj) - Mu) 

E(u) = -^ 7 ; 

u (l -^(u + iu) 



(29) 



It is of some interest to study Eq. I|29l) . with \i > 2 
and loT << 1. Using the imaginary Laplace transforms 
of [2£|, for u — > we get an analytical form for E(u), 
and with this analytical form plugged into Eq. I|27|). we 
recover Eq. l|25|l in the same limiting condition, namely 
the Laplace transform of ecos(a->t). This means that per- 
turbing the system in the infinitely aged condition, or 
keeping it under the perturbation influence, from the 
preparation to the infinitely aged condition, yields the 
same result. It is also interesting to remark that Eq. 
(|29|l establishes the effect of the step function pertur- 
bation in the case fi < 2. We must use Eq. (|29ll in 
the limiting case u> — > 0. This limit condition yields 
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E(u) = —(d^{u)/du)/(l — i>(u)), and with it, through 
the limit for u — > of i)(u), with /i < 2, and by means 
of Ea. l|27fl . the final result II(u) = e/u: The asymptotic 
limit is independent of the value of fi, and consequently 
of aging. 

The perturbation fit") = 0(t) cos(ojt), on the contrary, 
produces a response sensitive to the value of /i adopted. 
To show this important aspect let us make the assump- 
tion loT << 1, which, in the case fi > 2, turns Tl(u) into 
the Laplace transform of ecos(ojt), as shown by Eq. 1)25(1. 
In the case fi < 2, the assumption uoT << 1 allows us to 
obtain 



with Jn() denoting the Bessel function, and, more in 
general : 



n(t) 



cos (^/i — wi) 
T(^- l)(wt) 2 +M' 



(33) 



n(u) = ei?e 



(30) 



which means that the condition /i < 2 annihilates the 
response coherence in the long-time limit. 



Using the fractional derivative method of Ref. |2J we 
find Hi 



IT(t) = -eRe 



E 



fi - 1\ (-not) 71 1 t 
n J (n-1)! 



(31) 



-eRe [(fi - 1) F (2 - //, 2, lotf) e~ 



where -F (a, /?, z) is the confluent hypergeometric function 
|29|. In the special case /x = 3/2, of interest for the 
physics of blinking quantum dots 0, 0| , the method of 
Ref. yields Hg 



In conclusion, we derive the linear response of non- 
exponential renewal processes, without limitation to the 
stationary condition, and we show that aging, an un- 
avoidable consequence of /i < 2, has the interesting effect 
of annihilating the response to coherent perturbations, in 
the long-time limit. This result is obtained by means of 
the weak perturbation of a SME, which is equivalent to 
perturbing the corresponding CTRW. Along these lines 
it would be possible to consider also the case when the 
collision occurrence is significantly affected by perturba- 
tion. This would imply, however, a departure from the 
linear response. 



n(t) 



M ¥) cos (f 



(32) 
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